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Abstract — This paper describes the Probing In situ with 
Neutron and Gamma rays (PING) instrument, that can 
measure the subsurface elemental composition in situ for any 
rocky body in the solar system without the need for digging 
into the surface. PING consists of a Pulsed Neutron Generator 
(PNG), a gamma ray spectrometer and neutron detectors. 
Subsurface elements are stimulated by high-energy neutrons to 
emit gamma rays at characteristic energies. This paper will 
show how the detection of these gamma rays results in a 
measurement of elemental composition. Examples of the basalt 
to granite ratios for aluminum and silicon abundance are 
provided. 
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1. Introduction 

The Probing In situ with Neutrons and Gamma rays 
(PING) instrument [1] is a promising planetary science 
application of the active neutron-gamma ray technology that 
has been used successfully in oil field well logging and 
mineral exploration on Earth for decades. These techniques 
can be very powerful for non-invasive in situ measurements 
of the subsurface elemental composition on other planets. 
The objective of our active neutron-gamma ray technology 
program at NASA Goddard Space Flight Center 
(NASA/GSFC) is to bring instruments using this technology 
to the point where they can be flown on a variety of surface 
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lander or rover missions to the Moon, Mars, Venus, 
asteroids, comets and the satellites of the outer planets. 
PING differs from orbital gamma ray and neutron 
instruments in that it carries its own source of high energy 
neutrons: a high intensity Pulsed Neutron Generator (PNG) 
that produces fast (14 MeV) neutrons at a rate ~100x greater 
than neutrons generated via the interaction of galactic 
cosmic rays with planetary materials. This orders of 
magnitude higher neutron rate results in significantly shorter 
data acquisition times for the same precision. 

A critical measurement in the investigation of planetary 
bodies is the determination of the chemical composition of 
planetary crusts. In the course of planetary evolution, 
several processes modify the distribution of elements on the 
surface of a planet, namely mantle and crust differentiation, 
magma transport and its eruption onto a surface, and impact 
cratering. The abundance of certain elements associated 
with different types of geochemical behavior can provide 
important insight into the origin and evolution of these 
bodies. The global planetary composition is characteristic of 
the nature of its formation during the accretion from the 
solar nebula. In the course of planetary evolution, the 
elemental distribution is greatly modified by such processes 
as core and crust formation during differentiation and later 
magma formation and emplacement. A detailed discussion 
of the use of geochemical information in the study of the 
formation and evolution of planetary bodies can be found, 
for example, in [2,3]. 

PING will be an important tool for the above scientific 
investigations since it measures the bulk elemental 
composition as a function of depth without the need for 
drilling or digging into the surface. The ability to get this 
kind of information without drilling has many advantages: 
physically drilling into a planetary surface adds mass, 
power, volume and complexity to a mission and can pose 
great difficulties in fitting into constrained mission budgets. 


1 



Scientifically, the physical process of drilling into the 
surface causes problems with mixing the materials in 
various layers and can make depth measurements of 
elemental composition very difficult. PING can make 
subsurface measurements without disturbing the surface 
because of the penetrating nature of both the neutrons and 
gamma rays. 

PING is an adaptation of well-established techniques that 
have been used in the oil industry for decades. PING is 
capable of measuring in situ bulk elemental composition 
and density information over a ~ 1 meter radius and down to 
50 cm below a planets surface without the need for digging. 
While it is quite practical on Earth to take numerous 
samples of materials and send them to advanced laboratories 
for complete elemental analysis, it is extremely difficult to 
obtain samples from the subsurface of other planets in our 
solar system and bring them to Earth for similar analysis. 
The detailed analysis of the distribution of elements with 
depth provides useful information that can often answer 
questions about the chemical environment, the physical 
properties of the materials and the evolutionary processes 
that formed the geological structures. On Earth, the oil 
industry has led the way in the development of specialized 
instruments that can be used in the field to provide a depth- 
dependent elemental analysis and measurements of bulk 
properties that is almost as good as those performed in the 
laboratory. 

2. Active Gamma Ray/ Neutron 

INSTRUMENTS 

2.1 PING Operating Principles 

PING consists of three basic components: 1) a Pulsed 
Neutron Generator (PNG) that emits intense pulses of fast 
(14 MeV) neutrons to excite materials at and below the 
planetary surface, 2) a gamma-ray spectrometer to measure 
characteristic gamma rays from the excited elements, and 3) 
neutron detectors to measure the properties of the resulting 
lower-energy neutrons that reach the surface after scattering. 
As illustrated in Figure 1, when an extra-terrestrial surface 
is bombarded with fast neutrons from a neutron generator, 
the nuclei in the planetary material are excited and emit 
gamma radiation that is characteristic of isotopes of the 
elements present in the material. The intensity of these 
characteristic gamma rays can be directly converted to the 
abundance of the specific elements that 
emitted them, assuming an isotopic 
distribution typical of other bodies in the 
Solar System such as Earth. 

Additionally, fast neutrons incident on 
the planetary surface will lose energy 
with successive interactions within the 
material at a rate that depends on the 
material properties. By measuring the 
count rates of these slowed neutrons that 
reach the surface, we obtain information 
about density, hydrogen content, and 
subsurface layering. Active instruments 
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Figure 1: This sketch shows the PING neutron detectors, 
gamma-ray detector and PNG attached to the bottom of a 
planetary surface rover and illustrates the physical 
processes that produce the gamma rays and neutrons that 
PING detects. 

that stimulate planetary surfaces with fast neutrons and then 
detect both the gamma rays and neutrons emerging from the 
object’s surface thus provide valuable information about the 
subsurface material composition. 

These gamma rays are emitted at specific energies 
characteristic of each element via processes such as inelastic 
neutron scattering, thermal neutron capture, delayed 
activation, and natural radioactivity. Thus the elemental 
abundances in the stimulated material can be readily 
determined from analysis of the gamma rays emitted by 
these different processes. The probability of these processes 
occurring depends on the speed of the neutrons as they lose 
energy by scattering from fast (14 MeV) to epithermal (< 
500 keV) and thermal (< 0.4 eV) energies. Since the 
inelastic neutron scattering process requires the neutron to 
have a significantly higher energy than the reaction 
threshold in the nucleus (1-6 MeV), inelastic scattering 
events occur before the incident 14 MeV neutrons have had 
time to lose much energy. Thus the emitted gamma rays 
from inelastic scattering processes are produced earlier in 
time than gamma rays produced by other processes. 
Characteristic gamma rays resulting from thermal or 
epithermal neutron capture interactions appear a short while 
later. Following the inelastic scattering or neutron capture 
processes, the interacting nuclei may be left in a radioactive 
ground state that lives for some time (fraction of a second to 



Figure 2: Placement of timing windows relative to each PNG pulse. 
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many days) and then usually decays by beta decay, often 
accompanied by gamma ray emission from the daughter 
nucleus. Finally, there the naturally radioactive elements K, 
Th, and U that are common in solid bodies in the solar 
system and emit characteristic gamma rays from radioactive 
decay. No outside stimulation of these elements is needed 
for gamma-ray production. The intensity of the 
characteristic gamma-ray lines measured by the gamma-ray 
spectrometer can thus be used to infer absolute elemental 
abundances of the material. Elements accessible to PING 
include C, H, O, P, S, Si, Na, Ca, Ti, Fe, Al, Cl, Mg, Mn, 
and V. Excellent reviews of the physics of neutron - gamma 
ray techniques and how they are used for both remote and in 
situ geochemical analysis can be found in references [3-5]. 

2.2 PING Timing Advantages 

The fact that different nuclear processes occur at different 
relative times between fast neutron pulses can be used to 
significantly increase both the accuracy and precision of the 
inferred elemental composition of planetary materials. The 
use of a properly designed PNG provides great flexibility in 
setting coincidence gates to optimize both data collection 
efficiency and spectral resolution, resulting in data with 
increased precision per gamma ray detected. The diagram in 
Figure 2 and spectra in Figure 3 show how restricting 
gamma-ray data acquisition to specially timed coincidence 
windows between neutron pulses results in the separation of 
gamma-ray spectra by nuclear process. 

Bodnarik et al ., (2012) [6], explain in detail how the 
coincidence time window boundaries and the PNG pulse 
width and period can be optimized to separate gamma ray 
lines from these different processes. Our calculation in 
section 4 makes use of this ability to correct for spectral 
contamination from different processes as discussed when 
computing the rate of the 28 Si inelastic line at 1779 keV. In 


the case of the 1779 keV line there also happens to be 
delayed activation lines from both silicon and aluminum 
that result in a gamma ray at this same energy. Because the 
decay of these activated states are much longer than the 
PNG pulse, the 1779 keV gamma rays from activation will 
persist into the time of the following fast neutron pulse and 
will thus contaminate the inelastic spectrum at this energy. 
We can correct for this contamination because the activation 
line intensities are essentially constant throughout the PNG 
pulse period. The activation line strength can be measured 
during the activation window and then scaled by the 
window widths so that its contribution can be subtracted 
from the count rate in the inelastic window. This concept is 
the basis for generating the “corrected” data described in 
section 4. 

Figure 3 shows spectral data from a 6.3 hour run on the bare 
basalt in the experimental configuration described in section 
3. Here we see the spectra from inelastic scattering, neutron 
capture and activation well separated and showing different 
gamma ray lines and different background levels. 

3. PING Experiments 

3.1 Description of the GGAO Test Site 

We are testing the capabilities of our PING instrument 
prototypes at a unique outdoor gamma ray and neutron 
instrumentation testing facility located at Goddard’s 
Geophysical and Astronomical Observatory (GGAO) near 
the Goddard Space Flight Center main campus. A schematic 
aerial view of the test site is shown in Figure 4. This test 
facility allows us to operate PING on top of either of two 
large, assayed granite and basalt monuments, each 1.8 m x 
1.8 m x 0.9 m in size. Activation Laboratories Ltd. in 
Ancaster, Ontario, Canada, has independently measured the 
full elemental compositions of these Concord Gray Granite 
and Columbia River Basalt materials to the ppm level. The 



Figure 3. Spectra from different time windows: Gamma-ray spectra from a 6.33-hr acquisition using a HPGe detector on 
Columbia River basalt. 
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Figure 4. Aerial view of GGAO: This schematic of the 
outdoor gamma ray and neutron instrumentation testing 
facility shows the operations control building as well as 
the 46 m diameter safety perimeter surrounding the two 
existing 1.8 m x 1.8 m x 0.9 m granite and basalt 
monuments. 



Figure 5. PING experiment set-up: PING deployed for 
measurements on top of the basalt monument. The PNG is 
on the left, the HPGe detector is on the right, and 3 He 
detectors are between them. The data acquisition 
electronics are situated behind the basalt and are not 
visible in this photo. 

results of these assays are listed in Table 1. PING is 
remotely operated from a building 75 m from the 
monuments due to the radiation hazard from the PNG’s 14 
MeV neutrons. Underground power and communications 
lines connect the operations building to the test monuments. 


Further information about the test facility can be found in 
[7, 8]. 

3.2 The Importance of Using Both Granite and Basalt 

Both granite and basalt are igneous rocks containing large 
amounts of silicon and oxygen. While granite forms when 
magma cools slowly beneath a planet’s crust, basalt forms 
when deeper magma breaks through the crust, such as 
through a volcano, and then cools quickly on the surface. 
Because basalts are formed from deeper magma, they 
contain significantly more of the heavier elements such as 
iron and magnesium than granite. When PING is able to 
measure the elemental composition of a planet’s surface 
with sufficient precision to differentiate between igneous 
rock compositions, it will become a very effective tool for 
studying the composition of a planet’s mantle as well as 
local igneous processes. The first step toward this goal is to 
test PING with both basalt and granite and make sure that it 
can differentiate between these extremes in iron 
composition. As our techniques mature and we learn to 
optimize PING for sensitivity to higher iron content, the 
precision in PING’s material composition measurements 
will greatly improve. 

3.3 PING Experimental Setup 

Qualitative descriptions of PING results for various layering 
configurations such as with polyethylene and with the 
addition of materials such as titanium that are not found in 
either granite or basalt [1] have shown PING to be a 
powerful tool for subsurface elemental composition 
measurement. In this paper, we take the next step and begin 
the quantitative analysis of the elemental composition of the 
granite and basalt materials. We begin with a direct 
comparison between the bare basalt and bare granite 
configurations with PING resting directly on top of these 
large rock formations without any intervening layers of 
other materials. Figure 5 shows our experimental setup for 
the basalt. The identical configuration was used for the 
granite. Our PNG is a model MP320 from Thermo 
Scientific in Colorado Springs CO. Our gamma ray detector 
is an EG&G ORTEC Gamma-Gage coaxial n-type high 
purity germanium (HPGe) gamma ray spectrometer 
connected to a portable liquid nitrogen dewar as shown in 
Figure 5. This detector is used to measure gamma rays in 
the .05-10 MeV range with an energy resolution of 0.2% 
at the 1779 keV Si inelastic gamma ray line. In this 
experimental configuration, both PNG and HPGe detector 
are placed directly on top of the granite and basalt 
monuments so that they rest on the opposite corners, aligned 
along the diagonal as shown in Figure 5. The PNG fast 
neutron emission point is located 100.5 cm from the flat 
face of the HPGe detector and is 9.5 cm above the surface of 
the monument. The cylindrical axis of the HPGe 
spectrometer is located 14.0 cm above the surface. The 
relative locations of the PNG and HPGe were identical for 
the tests on the basalt and granite. Neutron detectors were 
placed between the PNG and HPGe gamma ray 
spectrometer as shown in Figure 5. 
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The PNG was set to emit an average of 0.4 xlO 8 neutrons 
per second with a pulse period of 1 ms and a 100 ps pulse 
duration. The gamma ray and neutron detectors were each 
read out in event-by-event mode by Lynx Digital Signal 
Analyzer units made by Canberra, Inc. using custom 
software. These data acquisition electronics were located 
below the surface of the formation so that the rock would 
provide shielding from the 14 MeV fast neutrons. The Lynx 
and PNG electronic clocks were synchronized to the PNG 
neutron pulses via a TTL output pulse generated by the 
PNG. This synchronization ensures 0.1 ps event time 
accuracy. 


granite densities are so similar (2.69 g/cc and 2.63 g/cc 
respectively). In fact, most of the remaining factors in 
Equation 1 also cancel if we take the ratio of the intensity of 
the identical energy line in the basalt and in the granite 
experiments. Since our ratio concerns the same gamma ray 
line coming from the same element A but in the basalt vs. 
the granite, the remaining factors cancel and we are left with 
the simple relation: 


^ ' A _ Basalt 
^ A -Granite 


f 
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R 


A _Basalt 


R 


A -Granite J 


( 2 ) 


The preliminary data presented here represent 4 hours live 
time on both the basalt and the granite. The measured dead 
times were both ~ 10%. 

4. Preliminary Results 

4.1 Comparing Data and Abundance Ratios 

As explained in section 2, the count rate for a gamma ray 
line that is due to a specific process on a single element is 
directly proportional to the abundance of that element in the 
material. Similarly, the equation below shows the linear 
relationship between the gamma ray line count rate ratio and 
the elemental abundance ratio. 


® A _Basalt _ 

R ^ 

^ A _Basalt 
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JD 
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Where 

(0 A_Basait C0 A_ G mni,e are the abundances of the same element 
A in the basalt and in the granite (in wt %) 

R A_Basalt R A_Gmni,e are the C0Unt rateS for a Specific gamma 

ray line from element A in the basalt and granite 
G A Basa i t is the interaction cross section for element A 

7a is the number of gamma rays with energy E A emitted 
per neutron interaction 

£(E a ) is the detector efficiency for gamma rays with 
energy E A 

Equation 1 is the result of the determination of the gamma 
ray count rate R for a given nuclear process involving a 
specific element. Equation 3 in Appendix A gives the count 
rate as a function of: 1) the macroscopic cross-section for 
both the bulk material and the specific element; 2) the 
number density of the specific element in the material; 3) 
the branching ratios for gamma ray emission of the specific 
line from a specific element; 4) detector efficiency, 5) PNG 
output, and 6) the absorption of emitted gamma rays as they 
travel through the rock to get to the gamma ray 
spectrometer. By taking the ratio of gamma ray rates for the 
same energy line from different formations, most of these 
factors cancel, especially considering that the basalt and 


The ratio of count rates from a specific gamma ray line from 
element A from basalt and granite is equal to the ratio of the 
abundance of element A in the basalt and granite. (See 
Appendix A for a complete derivation of Equation 2.) In the 
following sections, we compute these ratios for the 1779 
keV line from inelastic scattering on 28 Si and for the 2211 
keV line from inelastic scattering on aluminum 27 A1 and 
compare the ratio of rates to the ratios of the abundances of 
silicon to aluminum in the basalt and granite as given in 
Table 1. Table 2 contains the measured results needed for 
this comparison. 


Table 1: Granite and Basalt Assay Results 


Element 

Columbia River 
Basalt 

Concord Gray 
Granite 

Abundance (wt%) 

Abundance (wt%) 

O 

44.97 

48.59 

Si 

23.18 

34.59 

A1 

8.637 

7.32 

K 

0.955 

4.28 

Na 

2.448 

2.23 

Fe 

7.344 

1.00 

Ca 

6.618 

0.62 

Mg 

4.794 

0.17 

Ti 

0.898 

0.14 

P 

0.183 

0.09 

H 

0.087 

0.00 

Mn 

0.04 

0.125 

Rb 

0.035 

0.001 

C 

0.03 

0.03 

Ba 

0.026 

0.036 

Cl 

0.02 

0.04 

Zr 

0.014 

0.016 

Ce 

0.01 

0.006 

N 

0.01 

0.01 

Zn 

0.007 

0.01 

Sr 

0.006 

0.086 
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Table 2: Measured Gamma Ray Count Rates with PING 



Formation 

Egamma 

(keV) 

Net Count 
Rate (cts/s) 

Net Count 
Rate Error 
(cts/s) 

Basalt/Granite 

Ratio 

Basalt/Granite 
Ratio Error 

Independent 
Assay Ratio 

28 Si(n,n'y) 28 Si 

Basalt 

1779 

0.664 

0.037 

0.707 

0.101 

0.670 

28 Si(n,n'y) 28 Si 

Granite 

1779 

0.940 

0.123 

27 Al(n,n'y) 27 Al 

Basalt 

2211 

1.089 

0.019 

0.936 

0.023 

1.180 

27 Al(n,n'y) 27 Al 

Granite 

2211 

1.164 

0.021 


4.2 Ratio of Si in Basalt to Si in Granite 

It would be most convenient to simply divide the gamma 
ray line count rates to get the abundance ratios. However, 
when there are times where counts are produced by more 
than a single reaction on a single element, there is no longer 
a linear relationship between the counts in the peak and the 
concentration of a single element. To get the most accurate 
result for an element’s concentration it is necessary to 
remove the counts measured at times when they can be 
produced by multiple reactions or by different elements. 
One would like to have the 1779 keV peak that occurs 
during the high-energy neutron pulse be only due to silicon. 
However, there is also peak at the same energy that is due to 
the delayed activity of aluminum. Since delayed activity 
peaks are present at all times, if these counts were not 
subtracted from the peak measured during the high-energy 
neutron pulse, the derived elemental concentration would be 
much too high. 

Without the pulsing feature of our PNG and our use of time- 
resolved spectra, we would not be able to separate these 
processes. Because we can take gamma ray data during 
specific time windows between fast neutron pulses, we can 
isolate the contribution from the aluminum activation 
process and, after scaling by the ratio of time window 
widths, can subtract the contribution of this A1 activation 
line from the Si inelastic line count rate. Using a 90 ps 
inelastic window width and a 500 ps absorption/activation 
window width for the basalt data and 85 ps and 475 ps 
widths respectively for the granite data, we arrive at the 
following rates: for the 28 Si 1779 keV line, the gamma ray 
rate is 0.664 ± 0.037 cts/s from basalt, and 0.940 ± 0.123 
cts/s from the granite so that the count rate and hence the 
measured silicon concentration ratio is 0.707±0.101. This 
measured value is to be compared with 0.670 as determined 
by the elemental assay data shown in Table 1. We thus see 
good agreement between the measured and assay Si basalt 
to granite ratios; they agree to within one sigma. 

4.3 Ratio of Al in Basalt to Al in Granite 

To measure the aluminum abundance ratio between the 
basalt and granite, we will use the 27 Al(n,n’y) 27 Al inelastic 
scattering line at 221 1 keV. By fitting the aluminum lines in 
the basalt and granite spectra, we measure 1.089 ± 0.019 
cts/s for basalt and 1.164 ± 0.021 cts/s from granite, yielding 


a gamma ray count rate ratio, and thus aluminum 
concentration ratio of 0.936 ± 0.023. We note that this 
measured value is significantly lower than the 1.180 
aluminum ratio calculated from the assayed abundances 
values shown in Table 1. It is not surprising that this 
aluminum ratio agreement isn’t as good as the silicon 
because there is a fair amount of aluminum in the PNG and 
HPGe and neutron detectors that cannot be accounted for 
except by computer simulation. Thus, the aluminum 
measured by PING contains a constant value due to the 
PING equipment itself that must be added to the aluminum 
in the basalt or granite information and thus cannot be 
divided out in the ratio. Adding the same constant factor to 
both the numerator and denominator will reduce the ratio so 
that the result will be less than the expected ratio of the 
elemental assay data. This problem will be eliminated when 
our full Monte Carlo is running as described in Section 5 
below. 

5. Future work 

In the future, these data will be coupled with MCNPX [9] 
computer simulations to allow us to quantitatively 
determine the bulk elemental composition of the subsurface 
material for any solid body in the solar system. The 
MCNPX calculations allow an exact statistical calculation 
of both the energy and the time of a gamma-ray event 
detected in our gamma ray spectrometer. The calculations 
take into account all the details of neutron production and 
transport, all nuclear reactions on all elements present in the 
material, and all possible gamma rays that can be produced 
and their transport and detection at a specific point in space 
by a particular detector. Thus, the Monte Carlo calculations 
will provide a direct relationship between peak counts and 
elemental concentrations, limited only by the statistical 
count rate uncertainty and the systematic calculation 
uncertainties (generally less than 1%). 

Many of the parameters in the equations in section 4 and in 
Appendix A can only be determined using MCNPX 
computer simulations. The process by which gamma rays 
are generated and detected by the spectrometer is not 
analytically tractable. Indeed, not only are the neutrons 
emitted isotropically and can move around inside the 
formation material as they are moderated but they also emit 
gamma rays isotropically. A small fraction of the gamma 
rays will be going in the right direction to be intercepted by 
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the detector if they aren’t absorbed by the rock before they 
get there. Finally, many of the gamma rays that enter the 
detector will pass through it without being detected. The 
computer simulation traces all of these processes and 
includes the physics of neutron transport and interaction. 

Since landed planetary science instruments are limited by 
mass, power and volume, it is unlikely that it would be 
possible to use an HPGe gamma ray spectrometer for an 
actual mission due to the bulk and power required by its 
cryogenics system. Actual space-flight versions of PING 
would use spectrometers such as scintillators that have 
poorer energy resolution. We used the high resolution HPGe 
in these tests with known material samples so that we could 
best characterize the underlying spectra for the granite and 
basalt materials and thus better understand and characterize 
spectra taken with detectors with poorer energy resolution. 
The precision of elemental measurements is improved for 
increased exposure time. The chosen integration time for a 
measurement of a particular element is governed by the 
quantity of that element in the planetary soil and the 
measurement precision required. 

6. Conclusions 

We have measured the basalt to granite ratio for silicon with 
good agreement to the independent assay results. We have 
also measured the basalt to granite aluminum ratio with 
some success, especially considering that some of the 
measured aluminum comes from PING itself. The method 
described for comparing the elemental composition between 
the basalt and granite is not limited to the silicon and 
aluminum examples shown here. Using this procedure we 
can infer the basalt to granite ratio for the other major 
elements to arrive at a comparison of the total composition 
to within a few percent. The addition of Monte Carlo 
simulations will allow us to make independent quantitative 
measurements of elemental composition by performing the 
calculation of analytically intractable neutron and gamma 
ray transport factors that cancel in the ratio. 

This paper has shown PING to be a promising technique for 
in situ elemental abundance analysis of materials. Continued 
work to refine the technique and to build more flight-like 
prototypes will make it accessible to future flight 
opportunities for in situ measurements on any rocky body in 
the solar system. 
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APPENDIX A 

Derivation of the Ratio of Elemental Abundances as a Function of Count Rate Ratio 


We start with the equation describing the detector count rate from a specific gamma ray emitted by element A in any material. The 
following equation is appropriate for inelastic scattering interactions such as those that produced the silicon 1779 keV and 
aluminum 2211 keV lines used in the calculations shown in this paper: 

= 71 A® A NIyA A ) P y A 

where 

r a is the count rate for a specific gamma ray line from element A 
n A is the number density of element A in the material 
G A is the inelastic scattering cross section for element A 

A is the rate of neutron interactions and is proportional to the neutron generator output rate 

i 

Ja is the number of gamma rays with energy E A emitted per neutron interaction 

e(E ) 

' A ' is the detector efficiency at energy E A 

P 

Ja is the average gamma ray absorption probability for all gamma ray paths through the material 

Some of the above factors require Monte Carlo calculations due to complex geometry. However, they cancel out when you take 
the ratio: 

P \ _ Basalt ^ A _ Basalt ^ A ^^y A ^ ( P A ) Py A _ Basalt 

Pa _ Granite ^ A _ Granite ® A y A ^ ( P A ) PyA _ Granite 


Where n A 


P material ^ Avoga 


J A_ material r material y Avogadro 


Here A a is the atomic mass of element A, N Avoga d ro is Avogadro’ s number, p is the material density and 
CO A is the abundance of element A in each material. 

Substituting for n A we have: 


PA_Basalt ® A _ Basalt ^aP B asalt ^ Avogadro^ A, _Basalt^ yA^(P ' A ) PyA_ . Basalt 

PA_Granite \ ^ A _Granite j ^A P Granite^ Avogadro^ A _Granite^^ yA^(P ' A ) PyA _Granite 


Since the densities of basalt and granite are very similar and the probability of absorption P of the same gamma ray with energy E A 
will be the same in both materials, both of these factors cancel in the ratio. The input neutron factor N as well as Avogadro ’s 
number will also cancel. Solving for the abundance ratio we arrive at Equation 1 in the text: 


^ A _BasalJyA^P A ) 
® A _ Granite ^ y A ^ ( P A ) 


The detector efficiencies and branching ratio I are independent of the formation material and also cancel in the ratio. Finally the 
cross sections for gamma ray energy E A will cancel because we are considering the case where the gamma ray comes from the 
same element. 
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We are finally left with Equation 2: 


CO 


A_Basalt 


CO 


A _Granite 


R 

^A_ Basalt 

D 

V ^ A -Granite J 


Note that for gamma ray rates from lower energy neutron capture interactions, the formula is slightly different. The cross section 
for interaction must be represented as the fractional macroscopic neutron capture cross section in the material. We thus start with 
the following equation: 


f 


\ 


R 


A_ Basalt = 


n A -Basalt® A 


# elements 


^i -Basalt® i 


NI yA e{E A )P y 


A / 7 A _Basalt 


V i=i J 

Where all other factors are as defined on the previous page 
O ^ is now the neutron capture cross section for element A 


except that 
and 


# elements 

1 n, 


Basalt® i ls ^e total macroscopic cross section for neutron interactions 


;= l 

# e leme nts # e leme nts 

We make the excellent approximation that ^ Yl i Basait C) t ~ ^ Granite® ) so that this f actor cancels in the ratio. The 

i=l i = 1 

calculation then proceeds identically to that for the inelastic scattering interactions shown on the previous page. 
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